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ABSTRACT 

Color gradients in early-type galaxies contain valuable clues about their formation and evolutionary 
histories and mechanisms. We examine color gradients in 1,700 early-type galaxies in 159 galaxy 
clusters spanning a redshift range of 0.05 to 0.2. We find that color gradients strongly depend on 
the environment where galaxies reside, with steeper color gradients in poor rather than rich clusters. 
No dependence of color gradients on galaxy luminosity is found in either rich or poor clusters. The 
difference in color gradients can be explained by a change in the internal metallicity and/or an age 
gradient in these galaxies. Our results support a reasonable picture whereby young early-type galaxies 
form in a dissipative collapse process, and then undergo increased (either major or minor) merging 
activity in richer rather than in poor clusters. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: fundamental parame- 
ters 



1. INTRODUCTION 

Early-type galaxies (ETGs) have negative internal 
color gradients, with their stellar populations becoming 
bluer towards the galaxy outskirts. Different studies have 
shown that these gradi ents are mainly caused by a metal- 
licity variation (e.g. iPeletier . Valent iiin. and Jamesonl 
Il990|h As reviewed bv IPeletier et alJ (j!990t hereafter 
PDI90), color gradients can be used as an effective dis- 
criminant among galaxy formation scenarios. In a mono- 
lithic collapse picture, ETGs form by the rapid collapse 
of over-dense regions at high redshift. Galactic winds are 
expected to prod uce steeper metallicity gradien ts in more 
massive galaxies ljLarsonlll97l lCarlberdll984j) . On the 
other hand, merging is expecte d to produc e some dilu- 
tion of the gradients in galaxies ()Whiteil980l) , deleting or 
producing an inverse correlation among stellar popula- 
tion gradients and luminosity, and making the gradients 
less steep in denser environments. 

It is also unclear whether or not color gradients in 
ETGs depend on envir o nmen t and/or galaxy luminos- 
ity. iTamura and Ohtal (|2003l hereafter TaO03) found 
the mean color gradient of ETGs in a nearby cluster 
to be consistent with the val ue obtained for fie ld galax- 
ies by other authors (PDI90.IIdiart et aDl2002D. On the 
other hand, ITamura and Ohtal l)2000ft and lTamura et alJ 
(2000) studied color gradients in both cluster and field 
galaxies at intermediate redshifts, and found that the 
former seem to have less steep gradients. To date, no 
study of color gradients for galaxies in groups has been 
pe rformed. 

iKormendv and Diorgovsk] l)1989f) examined the de- 
pendence of color gradients on luminosity and other 
galaxy properties. While no clear trend with luminos- 
ity was found, they noted that the largest gradients oc- 
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cur in intermediate luminosity galaxies, and suggested 
that post-formation mergers on average diminish the pri- 
mordial gradients in the more luminous systems; this 
was also supported by their finding that galaxies with 
more anisotropic velocity distributions have weaker color 
gradients. TaO03 also showed that for cluster galaxies 
some evidence exists for a bimodal behavior, with very 
bright (L > L*) ETGs having gradients that steepen 
with lumi nosity, and fainter galax ies s howing the oppo- 
site tr end. lLa Barbera et ahl l)2004|) and lde Propris et alJ 
l)2004[) found that gradients in cluster ETGs do not 
change or can become more steep at fainter luminosities. 

In this Letter, we study the dependence of g — r color 
gradients in ETGs on galaxy luminosity and environ- 
ment. Color gradients are estimated by deriving struc- 
tural parameters, namely the effective radius, r e , and 
the Sersic index, n, for a sample of 1,700 ETGs in 159 
clusters 5 with different richnesses in the redshift range 
0.05 < z < 0.2. We assume a ACDM cosmology with 
n TO = 0.3, n A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. DATA 

The data are taken from the Palomar Abell Clus- 
ter Survey ijGal et alJ 12000^1 and were collected at the 
1.5m telescope at Palomar with a SITe 2048x2048, AR- 
coated CCD. The pixel scale at this detector is 0.368 
arcsec/pixel, yielding a 12.56' x 12.56' FOV. Data were 
ta ken in q, r, i Gun n— Thuan filters. We refer the reader 
to iGal et alJ pXXl for more details on data reduction 
and photometric quality. 

Galaxies used for this study were selected as fol- 
lows (see La Barbera et al. 2005, in prep, for details). We 
selected likely cluster members as those galaxies whose 
g — r colors are within ±a ■ <tcm around the color- 
magnitude relation of each cluster, where a cm is the 
observed dispersion of the red sequence. Field contami- 
nation was estimated using 34 blank fields observed with 
the same configuration and exposure times as the clus- 
ters. In order to keep field contamination as small as pos- 
sible without overly reducing the total number of ETGs, 

5 The list of clusters is available from the first author. 
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we use a = 0.5. From the 34 blank fields, we measure the 
mean number of field galaxies found in the same magni- 
tude and color range used to select galaxies in each of the 
clusters. Choosing a = 0.5, the fraction of field galaxies 
was typically ~ 10% for the entire cluster sample. We 
note that decreasing the value of a reduces dramatically 
the total number of ETGs without changing significantly 
the field contamination. For example, choosing a = 0.25 
decreases the number of ETGs by 40%, while the frac- 
tion of expected field galaxies decreases by only 1%. On 
the other hand, increasing the value of a makes field 
contamination significantly higher. For example, choos- 
ing a — 0.75, 1.0 or 1.5, the fraction of field galaxies 
increases to 16, 18 or 20% respectively. 

The value of <jcm was obtained for each cluster 
by fixing the red sequence slope to the value given 
bv iVisvanathan and Sandaeel l)1977|) at the correspond- 
ing redshift and by computing a 3<r clipped standard de- 
viation of residuals to the red sequence. We excluded 
galaxies with overlapping isophotes, which were defined 
using an isophotal level of 1.5<7 of the background stan- 
dard deviation. These objects could be interacting galax- 
ies at the cluster redshift or objects at different redshifts 
that overlap because of projection effects. This selection 
minimizes cases in which color gradients might be driven 
by interaction with nearby objects rather than by global 
cluster and/or galaxy properties. For each cluster, we 
also selected only galaxies brighter than a given magni- 
tude limit. We simulated galaxy images as a function 
of redshift and calculated the necessary S/N to recover 
the g — r col or gradients with system atic errors less than 
0.005 mag ( lLa Barbera et al.ll2002l hereafter LBM02). 
For each cluster, we set the magnitude limit such that 
this S/N is achieved in the radial range used (see below) . 
The typical value for the magnitude cut was r — 18.5. 

The above criteria yield a sample of 1,950 galaxies for 
which structural parameters were derived in the g and 
r bands. From these galaxies, we selected the 1,700 ob- 
jects w ith Sersic index n > 2.2, which are likely cluster 
ETGs (|Blanton et all 120031 For each cluster we mea- 
sure t he richness from DPOSS galaxy catalogs l|Gal et alJ 
|2003|) . Only 108 of the clusters (containing 1,135 ETGs) 
have well-calibrated DPOSS data. Additionally, a num- 
ber of systematic effects could introduce a bias into our 
results at higher redshift. These effects include the en- 
hanced detectability of rich clusters as opposed to poor 
ones, the faintness and smaller sizes of the galaxies in 
high-z clusters, and sampling of a different portion of the 
rest-frame galaxy spectra in the g and r filters. There- 
fore, for the richness dependent analyses we further re- 
strict the dataset to the 560 galaxies in 75 well-measured 
clusters with z < 0.15. 

3. DERIVATION OF COLOR GRADIENTS 

Surface photometry was derived by fitting galaxy 
stamps with PSF convolved Sersic models, as detailed in 
LBM02. For each stamp, a constant background value 
was fitted simultaneously with the model, achieving a 
typical accuracy of better than 0.1% in the background 
value. For each cluster and in each filter, the PSF model 
was constructed by fitting star images with a sum of 
three Moffat functions, giving a mean value of ~ 1 for 
the reduced x 2 of the PSF fits in the whole sample. 
Deviations of the PSF from a circular shape were also 
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Fig. 1. — For clusters with z < 0.15 and well measured richnesses, 
we plot the mean value of the g — r color gradients, V ff _ r versus 
the mean cluster richness, N ga i, (filled symbols). Each bin includes 
N = 140 galaxies. Empty circles are the mean values of V s _ r cor- 
rected for field contamination. Error bars denote la uncertainties. 
A small offset in richness has been added to separate the corrected 
and uncorrected data points. 



taken into account by expanding stellar isophotes in a 
cos and sin series. Details will be given in a forthcom- 
ing paper (La Barbera et al. 2005, in preparation). The 
best-fitting structural parameters (r G and n) in the g 
and r bands were used to derive the internal color pro- 
file for each galaxy, g — r(r), where r is the distance 
from the galaxy center. The color gradient, defined as 
V g _ r = S(g — r)/<5(logr), was estimated by the loga- 
rithmic slope of the color profile, performing a linear 
fit of g — r(r) vs. logr. For comparison with other 
works (e.g. PDI90), the fits were performed in the range 
r m iii < r < r max , with r min = 0.1 r c and r max = r , where 
r c is the effective radius in the r band. The typical 
value of r e for the galaxies analyzed here varies from 
~ 1.9"at z < 0.15 to - l.l"at higher redshifts. We 
also note that slightly varying (by ~ 30%) the values 
of r m ; n and/or r max changes the V g _ r values by less than 
0.01 mag, demonstrating the robustness of our results. 

4. RESULTS 

Galaxies were binned according to their parent cluster 
richnesses, each bin having the same number of galax- 
ies. As noted in §2, we use only the N = 560 galaxies 
in clusters at z < 0.15, which is the median redshift of 
the whole sample. Fig. ^ shows the mean color gradi- 
ent (filled symbols) as a function of the mean cluster 
richness N ga i. The figure clearly shows that V s _ r be- 
comes less steep for galaxies in richer clusters. The dif- 
ference in V g _ r between points with N ga i < 50 and 
N gal > 50 is -0.047 ± 0.008, and is significant at 5.9cr. 
In order to estimate the effect of field contamination on 
the V g _ r vs. Ng a i relation, we used a Monte-Carlo tech- 
nique. For each cluster, we measure N c galaxies which 
are likely to be ETGs. From the blank fields we then 
measure Nf galaxies which follow the same set of cri- 
teria for being in the cluster - these are the contam- 
inants. To account for possible magnitude-dependent 
differences between galaxy populations, the cluster and 
field galaxies are divided into three magnitude bins, and 
in each bin i we randomly select a total of N c j — Nf^ 
distinct galaxies which are expected to be cluster mem- 
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TABLE 1 

Mean color gradients for galaxies in 
clusters with different richnesses 
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TABLE 2 

Parameters of stellar population models for color 
gradients 



model 




ZUGyr) 


Zo/Zq 


T (Gyr) 


Vz 


V T 


Zi 


1.0 


11.4 


0.6 


11.4 


-0.22 





z 2 


1.0 


11.4 


0.4 


11.4 


-0.4 





TZ 


1.0 


11.4 


0.6 


8.0 


-0.22 


-0.15 


T 


1.0 


11.4 


1.0 


8.0 





-0.15 




redshift for high and low richness clusters (empty and filled circles, 
respectively). Error bars denote la uncertainties. Gradients from 
the literature for field and cluster ETGs are plotted as filled and 
empty symbols respectively, (see references in the figure). Curves 
represent the stellar population models for the color gradients. 



bers. This procedure was iterated 500 times for all of the 
clusters, re-computing at each iteration the mean color 
gradient in each richness range. Fig. plots the mean 
value of the corrected V 9 _ r values as the open symbols. 
The corresponding error bars are obtained by adding in 
quadrature the standard deviation of the background- 
corrected V g _ r values to the errors on the uncorrected 
gradients. The difference between the corrected points 
with N ga i < 50 and N gal > 50 is -0.047 ± 0.009, which 
is significant at ~ 5.5<r. We conclude that field contami- 
nation has no effect on our results. 



Fig. [5] shows the mean value of V 9 _ r in different 
redshift bins. The whole sample was divided into low 
{Ngai < 50) and high (N ga i > 50) richness clusters (here- 
after LRCs and HRCs), with similar numbers of galaxies 
in both subsamples. The gradients for each subsample 
are given in Table In Fig. [21 we also show optical color 
gradients from previous works for field (empty symbols) 
and cluster (filled symbols) ETGs at z < 0.2. The gra- 
dients were transform ed to V„_ r values using t he GIS- 
SEL03 spectral code ijBruzual fc CharlotJ 1200^) . as de- 
tailed in La Barbera et al. 2005 (in preparation). We 
see that our data are in good agreement with previous 
measurements of color gradients. Fig.[2]shows that color 
gradients slightly decrease with redshift, following two 
distinct trends depending on the cluster richness. To 
quantify these trends, we performed a linear fit of V g _ r 
vs. z for both LRCs and HRCs, using a least squares pro- 
cedure with V g _ r as the dependent variable. The uncer- 



tainties on the fitting coefficients were obtained by shift- 
ing points in Fig. |2 according to the corresponding error 
bars and repeating the fitting procedure. The slopes of 
the fitted lines are -0.52±0.I8 for LRCs and -0.32±0.2 
for HRCs. These values are negative at 2.9 and 1.6ct 
significance levels, respectively. We note that the dif- 
ference between V s _ r of LRCs and HRCs in Tab. H is 
—0.043 ± 0.008, which is significant at ~ 5er (cfr. above). 

To constrain the underlying age and metallicity gra- 
dients in the galaxies, we used four empirical stellar 
population models for V 5 _ r : two pure metallicity mod- 
els [Z\ and Z2), a pure age model (T), and a mixed 
age+metallicity model {TZ). Each model was con- 
structed by computing the difference in color between 
two simple stellar populations, describing the properties 
of stellar populations at the inner and outer galaxy radii, 
r m ; n and r max , used for the computation. Galaxy col- 
ors were obtained from the GISSEL03 code. The inner 
population was assumed to be old (T = 12 Gyr) with 
solar metallicity while the age and metallicity of the 
outer population were changed to reproduce the actual 
color gradients. In models Z\ and Z2 the outer pop- 
ulations are old, but their metallicities are changed in 
order to describe the color gradients in HRCs and LRCs, 
respectively. Model TZ was created to describe the gra- 
dients in poorer clusters by adding an age gradient to the 
metallicity gradient of model Z%, which reproduces the 
color gradient in richer clusters. In model T, the inner 
and outer stellar populations have the same metallicity, 
while the age gradient is chosen to reproduce the color 
gradient in rich clusters. Table reports the relevant 
parameters of each model, together with the correspond- 
ing metallicity and age gradients, V z = \og{Z Q /Zi) and 
V T = log(r o /r<) respectively. 

Fig. El plots the mean value of V 9 _ r for galaxies at 
z < 0.15 as a function of their absolute r magnitude, 
which was obtained by using k- corrections from the GIS- 
SEL03 code. Three samples are shown in the plot, cor- 
responding to galaxies in clusters with low, high, and all 
richnesses, respectively. We see that color gradients do 
not depend significantly on galaxy luminosity. 

5. DISCUSSION 

The dependence of color gradients on cluster richness 
can be explained by the fact that some physical mech- 
anism, related to galaxy environment, affects the age 
and metallicity gradients in ETGs. A pure metallic- 
ity gradient can account for the observed color gradi- 
ents, provided that Vz flattens from about —0.4 (model 

Z 2 ) in LRCs to 0.2 (model Z x ) in HRCs. These 

metallicity gradients can be compared with those esti- 
mated by previous works, which found them to be in the 
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tude, M r . Different symbols are used for HRCs, LRCs and for the 
whole sample, as shown in the lower left part of the figure. Error 
bars denote la uncertainties. 



ran ge [—0.3, —0-2], wi th a typical uncertainty of 0.1 (sec 
e.g. Ildiart et alJl2002j) . Both models Z\ and Z2, there- 
fore, are in good agreement with previous V2 estimates. 
The presence of overly shallow metallicity gradients in 
ETGs has often been invoked as a serious problem for 
monolithic scenarios of galaxy formation, sinc e these 
models predict values in the rang e [—1, —0.3] l|LarsorJ 
Il97l ICarlberelll98l IKawatalEfflll hereafter KAW01). 
Models Z\ and Z-i show, however, that color gradients in 
galaxies in LRCs are fully consistent with expectations 
from the monolithic formation scenario, but the consis- 
tency is only marginal for HRCs. A natural candidate 
for the variation of the mean color gradient as a function 
of environment is galaxy merging, which is expe cted to 
flatte n the stellar population gradient in galaxies \ White! 
1980). Our results support a scenario in which ETGs 
have undergone increased (either major or minor) merger 
activity in richer rather than in poor clusters. Models T 
and TZ, however, show that this interpretation is not 
unique. For example, the results for rich clusters can be 
equally well described by a metallicity gradient of —0.22 
or by an age gradient of —0.15 (model T). The ratio of 
these gradients is Vz/VV~3/2, w hich is simply the we ll 
known age- metallicity degeneracy l|Worthev et al. i f 19961. 
Model TZ shows that the richness dependence can also 
be explained by fixing the metallicity gradient of both 
HRCs and LRCs to that of model Z\, provided that 



galaxies in LRCs have a younger stellar population in 
the outskirts, with an age gradient of about —0.15. We 
note that a mild age gradient, of abo ut —0.1, would also 
be consistent with pre vious studies ijSaglia et alJ 120001 
ILa Barbera et all 12003). In a hierarchical merging pic- 
ture, galaxies tend to develop a disk by accreting gas 
from the surrounding halo, with this mechanism being 
inhib ited in high density environments ( Kauff mann et all 
1993). It would be conceivable, therefore, that a negative 
age gradient accounts for the steepening of the observed 
color gradients in LRCs. 

A troublesome issue for galaxy formation theories is 
the absence of any sharp correlation between color gra- 
dient and galaxy luminosity. In a monolithic formation 
picture, the same galactic wind mechanism which is re- 
quired to explain the color-magnitude relation of ETGs 
would produce a steepening of the metallicity gradient 
in more luminous galaxies (KAW01). Following KAW01 
(model B in their Tab. 3), we expect a V 9 _ r change of 
about —0.08 in the luminosity range of Fig. [21 which is 
in clear disagreement with our data. We note that a 
steepening of the color gradients as a function of galaxy 
luminosity is expected if ETGs form through the dissi - 
pative merging of disk systems (Bckki fc Shioval ll999). 
The absence of correlation between gradient and lumi- 
nosity could be explained by invoking (i) dissipationlcss 
merging, which could be more effective in flattening the 
metallicity gradient of massive galaxies, or (ii) a steepen- 
ing of age gradients in less massive galaxies, which could 
counteract the metallicity-luminosity relation. It is not 
obvious, however, which physical mechanism could drive 
this antagonistic behavior. The present data show that 
internal color gradients in ETGs probably result from 
a complex coordination of different physical processes. 
Crucial insights into these processes could be obtained by 
breaking the age-metallicity degeneracy, exploring color 
gradients for large samples of galaxies in the optical-NIR 
wavebands or at higher redshift. 
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